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Abstract 

Wavefunctions in chaotic and disordered quantum billiards are studied ex- 
perimentally using thin microwave cavities. The chaotic wavefunctions dis- 
play universal density distributions and density auto-correlations in agree- 
ment with expressions derived from a 0-D nonlinear cr-model of supersym- 
metry, which coincides with Random Matrix Theory. In contrast, disordered 
wavefunctions show deviations from this universal behavior due to Anderson 
localization. A systematic behavior of the distribution function is studied as 
a function of the localization length, and can be understood in the framework 

of a 1-D version of the nonlinear cr-model. 
PACS number(s): 05.45.+b, 03.65. G, 71.55.J 
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Complexity in Quantum Mechanics can arise from two sources : chaos and disorder. 
Non-Integrable systems which are classically chaotic are now known to be of relevance to 
a variety of atomic and nuclear phenomena [|TJ. Disorder is of fundamental importance 
in condensed matter physics, where it leads to the phenomenon of localization 0. For 
quantum chaotic systems, Random Matrix Theory (RMT) has been shown to provide a 
very good description of universal statistical properties of eigenvalue spectra |3||§]. Recent 
developments in the theory of disordered systems based upon nonlinear a-models using 
supersymmetry theory have led to the recognition that the extreme diffusive limit of 
disordered systems also behave similarly to quantum chaotic systems. These theories have 
made several quantitative predictions for correlations in mesoscopic systems, including the 
contribution due to localization. However these ideas remain to be tested experimentally. 
Experiments also have added importance as several systems of Quantum Chaos which have 
localization, may behave similarly ||, for e.g. the kicked rotator, which has been studied 
numerically shows deviation from RMT behavior 0. 

While the eigenvalue statistics of quantum chaotic systems has been experimentally 
tested, the eigenfunction statistics has not been done before. The principal reason is the 
lack of accessibility to wavef unctions. In atoms and nuclei, the nature of the wavefunctions 
only manifests indirectly in quantities such as transition rates. Experiments using microwave 
cavities which exploit the correspondence between the Maxwell and Schrodinger equations, 
are unique in that they allow direct measurement of the eigenfunctions, as well as eigenvalues, 
in model billiard geometries. Earlier such experiments have provided direct observation of 
[1(J, enabled precise tests of eigenvalue statistics and have demonstrated ability 



scars 



to study arbitrary geometries which are not accessible to numerical simulation. In this 
paper we use such experiments to study the influence of chaos and localization on quantum 
wavefunctions. 

The experiments were carried out using thin (d < 6 mm) cavities, whose cross-sections 
can be shaped in essentially arbitrary geometries. For these two dimensional cavities, the 
operational wave equation is (V 2 + k 2 )^> = 0. The details of the experimental method are 



described in ref. | 12| . Eigenf unctions were directly measured using a cavity perturbation 
technique which measures |\I/| 2 |T0|,|T2f . We study the manifestations of chaos and localiza- 



tion in terms of statistical properties of the eigenf unctions, such as the density probability 
distribution P(|\l/| 2 ) and the density auto-correlation P2(r) = < \^f(q)\ 2 \^>(r + q)\ 2 >, and 
the Inverse Participation Ratio (IPR), P2(r — > 0). 

The geometries studied were representative of integrable, chaotic and disordered systems. 
Among chaotic systems, besides the quarter Sinai billiard, we also study the Sinai-Stadium 
geometry introduced by us in ref. [jllj. The Sinai-Stadium has isolated periodic orbits (PO), 



and gives direct and exact agreement with RMT for the eigenvalue statistics, unlike the 
Sinai billiard and the stadium billiard for which the nonisolated PO contribution leads to 
deviations from universality (Tl| . Representative eigenfunctions for some of the geometries 
are shown in Fig.l. The 5.685 GHz (Fig. 1, top left) state of the Sinai-Stadium is similar to 
over almost a hundred other states of different energies that were observed, in that there are 
no visible scars that can be evidently associated with a PO. Indeed scarred states are few 
and the 7.370 GHz state shown in Fig. 1 (top right) is a rare coincidence with a whispering 
gallery PO. 

Localization effects were observed by fabricating billiards in which tiles were placed to 
act as hard scatterers. These geometries are obtained by placing 1cm square or circular 
tiles in a 44 cm x 21.8 cm rectangle at random locations (Fig. 1). (The locations were 
generated using a random number generator and the tiles placed manually). Earlier experi- 



ments [0] on the eigenvalue spectrum of similar disordered geometries had shown that these 
experiments exemplify textbook two dimensional electron systems without interactions to 
remarkable degree. To our knowledge, these geometries are difficult to study using numer- 
ical computation and hence the microwave experiments at present afford the only reliable 
means to study this type of disordered systems. More than one realization of the disordered 
geometry was experimentally studied, labeled Dl through D5 depending on the density of 
scatterers (Table 1). 

Sample eigenfunctions of D3, are shown in Fig. 1. Here again we have obtained nearly 100 



wavef unctions. It is evident that in these billiards, any association with classical structures 
such as periodic orbits is difficult to see, although of course these billiards are also chaotic 
and such association must exist in principle. Instead the most striking effect visible is 
localization, which is strongest in the 3.372 GHz state, but is also present in a weaker form 
in the 6.651 GHz state. The degree of localization can be varied either by changing the 
frequency window for a given geometry, or changing the density of tiles. Hence the mean 
free path / and also the localization length £ could be manipulated, the latter over almost 
two orders of magnitude. 

We first carry out an analysis of the chaotic wavef unctions. P(|\I/| 2 ) is shown in Fig. 2 
for the chaotic geometries - the Sinai-Stadium and the Sinai billiard. These are compared 
with the well-known Porter Thomas (P-T) distribution obtained from RMT [|i~4|j : 

P(|*| 2 ) = ^=ex P (-K) (i) 



The best agreement is shown by the Sinai-Stadium. The Sinai billiard displays slight 
deviations which we attribute to states influenced by bouncing ball orbits. In an integrable 
system, the probability distribution is often truncated, e.g. for a rectangle at \^\ 2 = 4, 
and is not universal. In contrast, the chaotic wavefunctions show a finite, although ex- 
ponentially vanishing probability of finding large intensities, in exact agreement with the 
P-T distribution, thus confirming their universal nature. Numerical simulations of random 
superpositions of plane waves was also done following ref. |1| and also shows very good 
agreement with P-T ||16|| . 



The density correlation P2(r) was also determined from the wavefunctions, and is shown 
in Fig.3. The angular brackets denote averaging over space Q , over angles between Q and 
r , and finally over several states. An important aspect of the present work is that our 
experimental results for the chaotic and weak disorder cavities (see later) are well fit by the 
functional form : 

P 2 (r) = l + cJl(kr) (2) 
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where k is the wave number. When computed for individual states, the result for P2(r) 
still has the functional form of Eq.(2), although the r — ► value is slightly different. However 
when averaged over several states, the resultant P2{r) is a robust quantity determined by 
the geometry alone. We further emphasize that we have checked and confirmed that i^O") is 
isotropic for the Sinai-Stadium. The asymptotic value of P2(r — > oo) = 1 can be understood 
since P2(r — > oo) — >< \^(Q)\ 2 >< \ty(q)\ 2 1 as the wavef unctions are completely 

decorrelated at larger distances. 

The data for the Sinai-Stadium yield c = 2 ± 0.1. The Sinai billiard yields c = 1.7 ± 0.1. 
An understanding of this value of c = 2 can be achieved by noting that in general, Pz(r — > 
0) =< |^(?)| 4 >= 1 + c. If x = |^(?)| 2 , then P 2 (r -> 0) = / °° x 2 P(x)dx. Hence from Eq. 
1 we get P 2 (r — > 0) = 3. Thus the value of c = 2 for the chaotic case are consistent with the 
P-T distribution in Fig.2. 

The functional form in Eq.(2) with c = 2 has been anticipated by Prigodin, et. al. PJI7| 
based upon a 0-D cr-model of supersymmetry. This theory is expected to apply to disordered 
systems in the diffusive limit, where the wavelength A << I (scattering length) << L, the 
cavity size. Results derived from this theory were earlier seen to coincide with RMT |Q. 
Therefore it gives P-T distribution for the eigenfunction components as in RMT. But it also 
gives coordinate correlations not available from RMT alone [16,171 . On the basis that for the 



diffusive limit, the randomization due to disorder of an infinite system is the same as that 
due to chaos at the boundary of a finite system, the results should apply to chaotic systems. 
Our experiments confirm both these predictions of the 0-D cr-model for P(|\I/| 2 ) and i^O") 
in the chaotic cavities. 

We now turn to an analysis of the eigenfunctions of the "disordered" geometries. We 
have divided the obtained wavefunctions into two regimes, (a) A > I, i.e. the wavelength 
is greater than the mean free path. Here one observes strongly localized regions (Fig.l, 
bottom left), (b) At higher frequencies, where A < I, the disordered eigenfunctions have 
the expected feature of having maxima apparently randomly distributed over the available 
domain (Fig.l, bottom right). In the experimental window from 2 to 8 GHz (15 to 3.2 cm), 



Dl is in regime (b), D2, D3 and D4 go from regime (a) to (b), and D5 is in regime (a). 

Fig. 4 displays the density probability _P(|\I/| 2 ) for D3 and D4, in regime (b) for data 
in the 6.5 to 7.75 GHz range. The disordered data are clearly not in agreement with the 
P-T distribution, displaying consistent and significant deviations at the high intensity end 
which are well beyond experimental error. These correspond to localized regions present in 
the disordered wavefunctions, in which the peak intensity is large compared to the average 
intensity. 

Recently, Mirlin and Fyodorov have been able to calculate the localization corrections in 
a quasi one dimensional wire using a 1-D nonlinear a- model of supersymmetry J7|. They have 
been also able to extend it to the full two dimensional case of the cavity in the delocalized 



regime [18]. For incipient localization, the results can be expressed as corrections to the P-T 



distribution given by: 



P(x) = /(x)exp(-x/2)/V27ra (3) 



where x = \^(Q) | 2 as before and f(x) = 1 + d(l/A — x/2 + x 2 /12) for small d. Now d 
depends on localization length £ and the density of disorder l/L through d = Lj (27c^)ln(L/l) 
for the two dimensional case ]T8l . This form is plotted in Fig. 3 for D3 and shows that the 



1-D cr-model gives an accurate description of the localization effects seen in the experiment. 
From the fit to data for Dl to D4 in the 6.5 to 7.75 GHz range, £/L are extracted and listed 
in Table 1. 

i^O") for D3 is shown in Fig.5 and is well described by the functional form Eq.2. The 
c value obtained from a fit is given in Table 1. This is consistent, within experimental 
error, with the values obtained directly from the density distribution in Eq. 3, by using 
IPR = J °° x 2 P(x)dx = 1 + c = 3 + d[L / (£ir)]lii(L / 1) . This higher value is therefore due to 
the high intensity peaks observed in Fig.l and the non-P-T distribution function in Fig. 4 
and depends on the localization length. 

In the strong localization regime (a) the deviations are even greater as the effects due to 
localization is even more pronounced. For strong localization, like in D5, stronger deviations 



in P(|\l/| 2 ) are seen which cannot be described by corrections to P-T form as in Eq. 3. 
However the data is described well by expressions for strong localization, for the wire case, 



in ref. 0. The form, P(x) = 8/(L/t)y(L/£)/2xK 1 (2^2x/(L/£)) fits well and a £/L value 
of 0.1 can be extracted for the data in the 3 to 6 GHz range, however some constants 
due to dimensionality considerations are expected to change. i^O") also does not obey the 
expression Eq. 2, as is evident from Fig.5. Instead spatial correlations die out faster, perhaps 
exponentially. Although the experimental data are available, further theoretical work needs 
to be done to compare the data in the strong and intermediate localization regime. 

An important perspective can be gained by examining the IPR for different eigenstates 
that were measured for both the chaotic and the D3 and D4 cavities, as shown in the inset to 
Fig.5. It is seen to fluctuate from one eigenstate to the other - such level to level variations 



are expected to have important characteristics [pl| which will be analyzed in future work. 
However the IPR for the chaotic cavity are confined to a band around the limiting value of 
3, with an average of 3.2. In contrast the IPR for the disordered cavity, show an increase 
towards the lower eigenstates, which tend towards the strong localization regime, showing 
the dependence of localization length on the frequency and the transition from regime (a) 
to (b). Of course deeper into the diffusion regime in which the localization length becomes 
effectively "infinite", at large / or small A, these also tend to 3, the value for the chaotic 
case. The dashed line represents the functional form for the IPR deduced above from eq.(3), 
and which extends up to the regime A < I, and which thus represents the small localization 
limit of the 1-D cr-model. 

In conclusion, we have been able to study the role of chaos and disorder in organizing 
quantum wavefunctions by means of electromagnetic experiments. The results demonstrate 
universal properties of chaotic wavefunctions, which also correspond to disordered systems 
in the extreme diffusive limit. They also demonstrate deviations from universality due to 
localization. The results appear to be consistent with a proposed notion of a (weaker) 
universality class |J for systems with localization, in which the statistical measures depend 
on one additional parameter - the localization length. The disordered billiards used here are 



the first experimental realization of such a system in which quantitative results have been 
obtained, and the dependence on the parameter £/L has been quantitatively demonstrated, 
enabling tests of the nonlinear a-model of supersymmetry theory, originally intended for 
mesoscopic systems. These results are of potential relevance to systems such as mesoscopic 
devices |HJ and acoustic and electromagnetic waves [20 . 
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FIG. 1. Representative eigenfunctions of the chaotic Sinai - Stadium billiard (top) and disor- 
dered geometry (bottom). 




FIG. 2. Density distribution Pd^j 2 ) for the Sinai-Stadium and the Sinai billiard compared 
with the Porter-Thomas form from RMT. While the Sinai-Stadium is in excellent agreement, the 
Sinai Billiard data shows slight deviations due to states influenced by bouncing ball orbits. 
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FIG. 3. Density correlation function P-2{r) for the chaotic geometries. The solid line represents 
Eq.(2) with c = 2 as expected from a supersymmetry approach. 




0.001 



0.0001 



5 10 . 15 20 25 

FIG. 4. Pd^j 2 ) for two realizations of the disordered geometries, showing deviations from 

Porter-Thomas due to localization. The dashed line is a fit to Eq. 3. The dashed-dotted line is fit 

to the strong-localization form in the text. 
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FIG. 5. P2(r) for the disordered wavefunctions showing oscillations similar to the chaotic data. 
The solid line represents Eq.(2) with c = 3.5 and kr rescaled by 0.96. (Inset) The inverse participa- 
tion ratio (IPR) vs. frequency for the D3 & D4 cavities (circles), and the chaotic cavity (squares). 
The dotted lines represent the limits of excursion of the chaotic data. The dashed curve represents 
the behavior expected on the basis of the 1-D a -model. 
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Cavity 


No of scatterers 


1, cm 


L, cm 




c 


Dl 


12 


9 


32 


10 


2.85 


D2 


27 


6 


31 


5 


3.14 


D3 


36 


5 


31 


4 


3.50 


D4 


36 


5 


31 


4 


3.54 


D5 


72 


3.5 


30 
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TABLES 



TABLE I. The various disordered cavities. 
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